j/as/?. ce-Af%z/? 


587 


NASA Contractor Report 159219 


NASA-CR- 1 592 1 9 

1180 Q0I0 { SS 

STRESS INTENSITY FACTORS IN TWO BONDED ELASTIC 
LAYERS CONTAINING CRACKS PERPENDICULAR TO AND 
ON THE INTERFACE - PART II. SOLUTION AND RESULTS 


Ming-Che Lu and F. Erdogan 


LEHIGH UNIVERSITY 
Bethlehem, Pennsylvania 18015 


NASA Grant NGR 39-007-011 
January 1980 


NASA 

National Aeronautics and 
Space Administration 

Langley Research Center 

Hampton, Virginia 23665 



FEB 2 1 1980 


BANGLEY RESEARCH CENTER 
LIBRARY, NASA 

VIRQIN'* 





STRESS INTENSITY FACTORS IN TWO BONDEO 
ELASTIC LAYERS CONTAINING CRACKS PERPENDICULAR 
TO AND ON THE INTERFACE - PART II. SOLUTION AND RESULTS 

by 

Ming-Che Lu and F. Erdogan 
Lehigh University, Bethlehem, PA 18015 

Abstract 

The analysis of the title problem was given in 
Part I. In this part the numerical method for solving 
the problem is described and the stress intensity factors 
obtained from the solution for various crack geometries 
are presented. 

1. Introduction 

The formulation and analysis of the plane problem 
for two bonded infinite dissimilar elastic strips which 
contain cracks of various configurations is presented in 
Part I of this paper [1], The problem is intended to approx- 
imate a composite beam or a plate having cracks perpendi- 
cular to and on the interface of the two layers. The 
crack geometries and some of their limiting cases which 
have been considered in the analysis are shown in Figure 1, 
This part of the paper is devoted to the solution of the 
problem for various typical crack geometries and to the 
presentation and discussion of the results. (*) 


(*) This work was supported by NSF under the Grant ENG 
78-09737 and by NASA-Langley under the Grant NGR 39- 
007-011 . 
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2. The Sol ution 

In Part I the problem was shown to reduce to a system 
of four singular integral equations for the unknown func- 
tions <j>.j , . . ,$4 which are defined by 


■jl- v.(x.,+0) - *,(*,) , (1*1,2) (!) 

^■[v^+O.y) - v 2 (2h 2 -0,y)] * i> 3 (y ) , (2) 


«[Ui(+0,y) - u 2 (2h 2 -0,y)] = ; 4 (y! , (3) 


where and v . » (i = l , 2 ) are respectively the x and the y 
component of the displacement vector in the strips 1 and 
2, and the general notation is shown fn Figure 2, The 
system of singular integral equations are given by (19)- 
(21 ) of Part I, The integral equations are solved by using 
a Gauss-Chebyshev or a Gauss-Jacobi integration formula 
whenever possible* and a combination of a Gauss-Jacobi 
integration formula and the method of Jacobi series if the 
simple Gaussian integration methods are not applicable* 
in order to apply these techniques it is necessary to 
normal ize the supports of the integral equations (19)-(21 ) , 
These equations are of the following general form (see 
Figure 2 ) : 




4 f bj ’ 

i C | h i j ( y ’S)<l>j(s) ds + y s j j i>j(y)l = p i (y ) , 

i = 3 , 4 , a i <y<b • (b) 

where a 4 - 'a,, b 4 = b 3 * the constant y is given by ( 22 ) 
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of Part I, P i ,..,P 4 refer to the crack surface tractions 
and are material constants (see equations 9-11 

of Part I). In the case of imbedded cracks the kernels 
h vl , h 22 , h 34 » h 43 have a Cauchy-type singularity and the 
remaining kernels are bounded. For cracks intersecting 
the boundaries or each other all related kernels have 
generalized Cauchy-type singularities. 

For the purpose of numerical solution the following 


normalized quantities are introduced: 


r = 


2s 

b .+a • 
3 J 

b r a j 

b .-a . 
j J 

2x i 

b . +a . 

b . -a • 

l i 

b i- a i 

2y _ 

b 3 +a 3 

b 3" a 3 

b 3" a 3 


( 6 ) 


, (a.<x.j<b. , -l<r<l, i *1,2), (7) 


( 8 ) 


f.(t) = <f-i (S ) . (9) 

J J 

h^jCx.jS) 58 H^j(r»t')» (a.<x-<b^, a j <s<b j’ ~l < ( r > t ) < l> 

i =1 ,2 ; j = l . ,4) , (10) 

h.jj(y»s) = H.j(r,t), (a 3 < y < ^ ) 3 » a 4 =a 3’ b 4 =b 3> 

a.<s<b., -l<(r,t)<l, i=3,4; j=l,..,4), 

J J 

(ID 

p i (x i ) = P i (r), i = l ,2, a.<x.<b., -l<r<l, (12) 

P|(y) = P,-(r) , i = 3 , 4 , a i <y<b i , -l<r<l. (13) 
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The system of integral equations (4) and (5) may then he 
expressed as 


4 


E 

3 = 1 


b .-a . 
_J 1 


4 

£ C 

3 = 1 



1 

| H • j (r,t)f j (t)dt = M^P|(r) , f=l ,2 , -l<r<l , 

-1 (14) 

1 

’ H 1j (r,t) f j (t) dt + yS^fjCr)] = M.P.fr), 

- 1 

i = 3,4 , -1 <r<l . (15) 


Also, defining 

f(r) * f 4 (r)+i f 3 (r), P(r) = P 4 (r)-i P 3 (r), Cl6a, b) 

and combining the two Integral equations given by (15) 
we obtain Csee equations (20), (21), (26) and C27) of 
Part I) 

1 1 2 

£ f W dt - V Hr! + I Kj(.r»t)fj(t)dt 

-1 -1 

1 

+ f [K 3 (r,t)f(t) + K 4 (r,t)f(t)]dt = M 3 P(r), 

-1 

-1 <r<l , (17) 

where the kernels K-,,..,K 4 are complex functions. 

Referring now to equations (24) and (28) of Part I 
and (9), (6) and (16a), the solution of the integral 
equations (14) and (17) may be expressed as follows: 

F .(t) 

f . (t) = — av — r. » 0<Re (a i ,0.. )<1 , 3 = 1 ,2 , (18) 

J (1+t) j (l-t) 3 J J 


f(t) = ^ TZ ’ 0<Re(a 3 ,6 3 )<l . (19) 

(1+t) 3 (l-t)^ 3 
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where for some typical crack geometries the characteristic 
equations to determine the constants ctj and 3 j > (3*1 >2,3) 
are given in Part I. For example, for the imbedded cracks 
we have 

a^j - 3-j ~ a 2 ~ s ^ » a 3 = \ " 1a) » ^3 = ^ + iw> 


“ = H 109 ( W> ' 


( 20 ) 


Note that the integral equations (14) are of the first 
kind and the constants and 3j > (j = l,2) are always real. 
Consequently, the related Chebyshev polynomials T n ( t ) (cor- 
responding to aj = 0.5 = 3j) or Jacobi polynomials 
P (-^j»- a j)(t) have n real roots t k , -l<t k <l. Therefore, 
in this case a Gaussian integration formula can be developed 
to evaluate the singular integrals and the integral equations 
may be replaced by a system of algebraic equations for the 
unknowns F^tjJ, j-1.2, T n (t k )=0, (or pJ' 6 r‘V(t k ) = 0), 
k=l,..,n (see, for example [2]). On the other hand (17) 
is a singular integral equation of the second kind having 
a dominant part with complex coefficients. Therefore, the 
constants a 3 and 3 3 will always be complex and the roots of 
the related Jacobi polynomials P n ^“ e j’“ a 3 (t) will have 
complex roots. Since these roots are not on the line 
of integration - 1 < t <1 , it is not possible to develop a 
Gaussian integration formula to evaluate the dominant part 
of the integral equations in terms of a discrete set of 
unknowns F(t k ), k=l . ,n. A convenient method to solve 
this integral equation is the method of related orthogonal 
polynomials described in . [3] (see also [2] and [4] for 
applications). Referring to (19) and observing that 


w ( t ) = (l-tro+tr , (a = - 3 3 * 


3 ■ - a 3 ) >. 


( 21 ) 


is the weight function of the Jacobi polynomials p|> a *^(t), 





In order to solve the integral equation the unknown .unc- 
tion f(t) is expressed as 

f(t) = Z c.P. Ca » e) Ct)w(t) 5 ( 22 ) 

o 3 J 

where c 0 , . . ,c N are unknown complex coefficients. Using 
now the following property of the Jacobi polynomials 


1 


X f P. (a>6) (t)w(t)^r - Y P. {a,3) (r)w(r) 

7T1 3 o-r J 

2T(1-Y 2 ) 1/2 <»-> * - 1<r<1 - 

the integral equation (17) may be written as 

1 

1/2 „ D (-a,-6) 


(23) 


z i(i-v 2 ) ,/z c 1 p;:?-* ,,, (r) * j 

j = 0 41 J J -1 


Z K,(r,t)f,(t)dt 
«| J J 


" f 1 [c.K,(r,t)p'“ lB) (t)w(t)+c J K 4 (r,t)P< c ‘’ B) (t)i(t)]dt 

3=0 J -T ° 


= M 3 P(r) , -1 <r<l . 


(24) 


Using now the following Gaussian integration formula 
which corresponds to the weight functions defined by (18) 
[5] 

9^}— q— = Z W k g(t k ) , ( 25 ) 

after evaluating the integrals (24) may be expressed as 


(-a,- 8) 


2 n 


j = 0 


J r( 1-Y*J 1 'Vi-r (r) + Z - K m (r ’ t k )Vl k F m (t k ) 


2i 


m-.l k-1 


+ Z [ c j Q 3 j ( r ) + CjQ 4 j(r)] 3 M 3 P(r), -l<r<l. 

j 0 


(26) 
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Note that 


w(t) = (i + t)’ a3 (i-t)’ e3 = (|^|) 1w n-t 2 ) 


- 1/2 


(27) 


Hence the functions Q 3 j and Q 4j . may be evaluated quite 
simply by using the following Gaussian integration formu- 
la [5] 

.... .. TT _ / j. \ . _ ...... ?! 

2n 

(28) 


9 .-V 2 n 7 r , x . ir(2k-l ) 

p(t)(l-t 2 ) dt = I : J p(t k ), t k = cos 


■1 


In order to reduce (26) to a system of algebraic 
equations for the unknown constants Cj and F m (t k ) one may 
use a simple weighted residual method. In this case, it 
is clear that the related orthogonal polynomials are the 
Jacobi polynomials a * 3 ^ ( t } , n=0,l,. .. . Thus, using 
the orthogonality relation 

f l (a, 6) (a,$) 

’n 


P_ (t)P_ ' ’ (t)w( t)dt = { 


0 , n f m 
6 n (a,6) , n=m , 


(29) 


■ , 0 \ 2 a+3+1 r(n+a+1 )T ( n + B+1 ) 

e n^ a ’ 3 ^ " (2n+a+B+l ) (n+a+B+1 )n: 


multiplying both sides of (26) by 
-a -6 (-a, -3) 

( 1 - r ) (1 +r) (r) , ^ = 0,1,.. ,N-1 

integrating in (-1,1) we find 


1/2 2 n 

ir ( l - Y 2 ) c e+1 e z (- a ,-e) a imk F m (t k ) 


2i 


+ l C b £j c j +<1 > (£=0,1 . , N — 1 ) . 

J 0 


( 30 ) 


-?- 


and 



Again, the constants a, b, and d may be evaluated by 
using the Gaussian integration formula (28), 

Similarly, the integral equations (14) may easily 
be reduced to a system of algebraic equations in F m ( t ^ ) 
and ty The first two terms in (14) may be expressed 
in terms of- F (t fe ) by using the integration formula (25) 
and (18) [2] giving 


' H fj (r s ,t)fj(t)dt S " w k H tjt r 's*t k ) F j (1: k ) • 

-1 K 

i = 1 ,2 , j*l ,2, s - 1 , . . ,n-l (31 ) 


Where t k and r $ are the roots of the 


corresponding 


Chebyshev or Jacobi polynomials. The last two terms on 
the left hand side of (14) Which contains f 3 and f 4 
involve real and imaginary parts of the integrals 


f H ij (r s’ t)P / a,3)(t)w(t)dt ’ 

-1 

which may easily be evaluated by using the integration 
formula (28) . 

It should be observed that by using the procedure 
outlined above the integral equations (14) and (17) have 
been reduced to 2n-2+2N real algebraic equations in 
2n+2N42 real unknowns F m (t k ) , (m=l ,2 ; k*l , . . ,n ) and the 
real and imaginary parts of c ^ , ( j=0 , , . ,N) , The remain- 
ing four equations are provided by four additional condi- 
tions corresponding to the single-valued ness of displace- 
ments, the requirement of boundedness of the displace- 
ment derivatives, and the relations between the displace- 
ment derivatives for intersecting cracks (see Part I ) . 

For example, for an imbedded crack 0<a.j <b.| <2h-j (see Fig- 
ure 2 ) 

r b] n 

(s)ds=0, or E | F 1 (t k ) = 0 
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(ir\ 

and for an edge crack 0<a^<b^ = 2h-| , F-j(l) = 0, etc. 1 
The details of the numerical solution for each typical 
crack geometry may be found in [6], 

From the definitions given by (1), (2), (3), and (9) 
it is seen that once the functions f-,,.. ,f 4 are determined, 
one may obtain the desired crack surface displacements 
by means of routine integrations. 

3. Stress Intensity Factors 

The stress intensity factors were defined in Section 
5 of Part I. In Part I it was also shown that the stress 
intensity factors may be calculated directly from the 
asymptotic expansion of the displacement derivatives 

*4*4 around the crack tips. Thus, aside from a multi- 
plying factor, the constants Fj(+1), (j = l>,.»4) give the 
stress intensity factors. For example, from equation 
(62a) of Part I and (6), (9) and (18) it follows that 
for an imbedded crack 0<a 2 <b 2 <2h 2 the stress intensity 
factors may be expressed as 

k C a 2 ) * TT7Z f 2 ( - 1! •'U> 2 -a 2 77? . 

k(b 2 ) = - ^ F 2 (l) /(b 2 -a 2 )/2 . (33a, b) 

Similar expressions can be obtained for all other crack 
geometries. 

4. Results and Discussion 

The particular crack configurations for which the 
problem is solved numerically are shown in Figure 1, 

The results are given in Figures 3-36. With the exception 

(*) Tn this case f -| C+l ) has the indefinite form 0/0 and , 
if needed, may be evaluated as a limit, 
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of a few crack geometries for which crack surface displace- 
ments are calculated, the results given in this paper refer 
to the stress intensity factors. Generally the results 
presented in the figures are self-explanatory, Hence a 
detailed discussion of each figure dees not seem to be 
necessary. 

The results given in this paper are obtained for self- 
equilibrating crack surface tractions . If the externa] 
loads are applied to the layered material at locations suf- 
ficiently far from the region of cracks, then the crack 
surface tractions in the perturbation problem would be 
uniform* For example, if the medium is loaded in tension 
parallel to the y-axis away from the crack region the 
crack surface tractions are constant and are related by 


(1-v 
' £ 


Mil 

1 


( 1 -V2 ) P2 


?3 



for plane strain and 


Pi . El 

E 1 E 2 


P3 



(.35) 


for plane stress . In this paper pply plane strain case 
is considered. 

Figure 3 shows the stress i nteps i ty factors for a 
simple imbedded crack of length f|.g, '? hg? Note that as 
the crack tip b 2 approaches the interface (i .e. , as 
k(b 2 ) tends to zero for u 2 <p^ and tP infinity for I4g >v l * 
This wel 1 -known behavior i s due to the f§et that for b 2 =2h 2 
the power of the stress singularity $2 1 § greater than 
0.5 if u 2 > h 1 anc * less than 0.5 if u^v-j [see> for example 
[ 7 ] ]. For this case the definition of the stress inten- 
sity factor and the behavior of the stress state around 
the crack tip are given by the equation? (65) and (57) of 
Part I. For the material combinations psed in Figure 3 
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0 2 =0.624348 for y^y-, and & 2 = 0.385339 for y 2 <y r 
Also note that as the crack tip a 2 approaches the free 
boundary as expected, k(a 2 ) tends to infinity, 

Stress intensity factors for an edge crack are shown 
in Figures 4 and 5. In this case too note that as the 
crack tip b 2 approaches the interface k(b 2 ) tends to 
zero for u 2 <y-| and to infinity for v 2 >y'j . Also note that 
as the crack length decreases the stress intensity factor 
ratio approaches 1.586 which is the value obtained for 
the semi-inf 1 ni te plane having an edge crack of length 2^ 2 
Figure 6 shows the effect of stiffness ratio for a fixed 
crack geometry. Note that as y-j/y 2 ■+■ 0 the stress inten- 
sity factors approach the value for an edge-cracked strip 
under uniform tension k(b 2 ) = 3,99 p 2 vXJ\ The effect of 
thickness ratio is shown in Figures 7 and 8, The results 
for the crack terminating at the interface are shown in 
Figures 9-11 as well as in Figures 7 and 8 (the curves 
corresponding to -^/h 2 = l). 

Figures 10 and 11 show the stress intensity factors 
in a layered plate containing an edge crack and subjected 
to uniform bending away from the crack region. In the 
uncracked composite plate the relevant stress is given by 
(see Figure 2 and insert in Figure 10} 

a 2yy (x 2 ,0) = -p 2 (x 2 ) = -p 2 (l-x 2 /c 2 ) (36) 

c 2 = CE 1 (h2 + 2h 1 h 2 )+E 2 h2]/(E 1 h 1 +E 2 h 2 ) (37) 

where x 2 =c 2 determines the location of the neutral axis 
and the constant p 2 is the magnitude of the stress at the 
surface which is related to the bending moment M by 


Vo = 


3c 2 M 


2(2h ] -c 1 ) +3c^j -c 2 


2 j » c i — 2 ( h i h 2 ) - c 2 


( 38 ) 
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One may note that since the powers of stress singularity 
are different, the direct comparison of the stress Inten- 
sity factors for the two cases shown in Figure 11 would 
be meani ngless . 

The stress intensity factors in bonded layers con- 
taining two (.col 1 inear ) edge cracks and subjected to ten- 
sile loads away from the crack region are shown in Figures 
12 and 13. Note that in this case the stress intensity 
factor in the layer with the smaller stiffness is generally 
smaller than that in the layer with greater stiffness, 
decreases with increasing crack length, and eventually 
becomes negative. This, of course, is due to the "bending" 
of the notched composite layer. In the figures the nega- 
tive stress intensity factor is shown by dashed lines. 

If there are no other external loads to offset this effect 
through superposition, the results given in the figures 
for k<0 are not valid. 

To give some idea about the effect of Poisson's ratios 
on the stress intensity factor in the case of a broken 
layer Figure 14 shows the results for varying v-j and fixed 
values of v 2 "0.3' and v-j*2v 2 * lt is important to note that 
in this case the power of stress singularity £ 2 =a is also 
a function of v.,, which is given in Table 1. Thus, even 
though the figure shows the stress intensity factor slightly 
increasing with v-j, since a is a decreasing function of 
v i , the intensity of the stress state at the crack tip 
would actually decrease with v-j (see equations (65) and 
(67) of Part I). 

Table 1. Power of stress singularity 6 2 = a for a crack 
terminating at the interface, >>2*0.3, y-j *2y 2 . 

1 0 0.05 0.10 0.15 0,20 0.25 

a 0.483251 0.476747 0.469679 0.461967 0.453516 0.444205 

1 0.30 0.35 0.40 0.45 0.50 

a 0.433891 0.422389 0.409466 0,394818 0,378041 


- 12 - 



Stress intensity factors for a crack crossing the 

interface are given in Figures 15-19. The results for 

an imbedded crack are shown in Figures 15 and 16, Note 

that as the crack tip a 2 approaches the interface, due 

to again the change in the power of stress singularity, 

the stress intensity factors become unbounded. Also as 

a 2 -*-0, k(a 2 )-H». The results for the edge crack (i.e., for 

a 2 =0) are given in Figures 17-19. As h^0,2 h 2 , the 

crack tip b^ approaches the free boundary and consequently 

k ( b 1 )-»•«. Figures 18 and 19 show the stress intensity 

factors k and k at the interface. It should be noted 
x xy # 

that even though the magnitudes of k x and k xy increase 
rapidly with decreasing net ligament thickness 2 (h.j -0 . 2h 2 ) , 
h 1^2 =0,2 1S not an asymptote * x and k X y' Figure 20 
shows some sample results for the crack surface displace- 
ment . 

An important special case of the crack crossing the 
interface is the stress-free end problem for two bonded 
semi-infinite plates or beams. In this case the problem 
may be solved by letting a.| = 0=a 2 , b -j =2h ^ and b 2 =2h 2 (see 
Figure 2) and properly treating the singular behavior of 
the solution at the irregular points x 2 =0, x 2 =2h 2 , x^-0, 
and x -j = 2 h -j (see Part I). The problem may be one of ther- 
mal stress, residual stress, or mechanical loading. The 
technique for isolating the perturbation problem in which 
self-equilibrating end tractions are the only external 
loads is relatively straightforward. Invariably one 
imagines an infinite plate (- oo <y< co , -<»<z<» , 0<x-|<2h^, 
0<x 2 <2h 2 ) under a given set of applied loads and calculates 
the stresses at y=0 plane. The needed tractions in the 
perturbation problem are the equal and opposite of these 
stresses. For example if a clamped plate is pulled in 
z-direction or heated uniformly, the perturbation problem 
would have to be solved under a system of self-equilibrating 
tractions on the boundary y=0 which consists of constant 
normal stresses p-, and p 2 and the linearly distributed 
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stress coming from a bending moment M which are related 

by 

p ^ h i + ® M 55 2 p -j h-j (h-j+hg) . (39) 


The stresses resulting from M would be obtained from 
equations such as (36)-(38). Figure 2i shows the results 
of such art example. 

Ih the stress-free end problem* after solving the 
integral equations and obtaining the density functions 
<j, 1 and <J >2 the relative displacements in y-direction on the 
boundary y=0 may easily be evaluated (See equation (1 ) ) 


from 

V^x^O) = v Q 
V 2 (x 2 »0) = v Q 


X 1 

» 

4>1 ( s ) d s , 

D 

2h 2 

» ■" 

<h 2 ( s ) ds * 

x 2 


(40a ,b) 


Where V Q * v ] (o,o) = v 2 (o,o) represents a rigid body 
motion. Figure 22 shows some sample results for the dis- 
placements . 

The res ul ts for a uniformly pressurized T-shaped 
crack are shown in Figures 23-25. Note that due to pres* 
sure p ^ for Small values of the crack length 2£ 2 the 
stress intensity factor k (a 2 ) becomes negative (Figure 24). 

Figures 26-29 show the resul ts for bonded 1 ayers 
containing a T-shaped crack and loaded in y-directioh away 
from the Crack region for which Pg^P^'O 4hd P 2 is constant. 

The results for a symmetric cross-shaped crack are 
given in Figures 30 and 31. The external load in this 
example is the tension in y-directiort away from the crack 
region. Corresponding results for a 2 =0* i.e>, for a broken 
layer are shown in Figures 32-34. The Crack surface 
pressures Used in the example given id Figures 32 and 33 
correspond to uniform tension ih x direction as well as 
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tension in y direction away from the crack region. The 
results for tension in y direction only are given in 
Fi gure 34 . 

The results for an edge crack in a bonded semi- 
infinite plate are shown in Figures 35 and 36. The ex- 
termal load in these examples correspond to a uniform 
tension in z- direct ion or uniform heating or cooling of 
a clamped plate. Note that for this crack configuration 
the dominant stress intensity factor is the shear compo- 
nent k- 2 • From Figure 36 one may also observe that for 
bonded plates with equal thicknesses the magnitudes of 
the stress intensity factors are rather small. 

Finally, it should be pointed out that the accuracy 
of the numerical results given in this paper is not uni- 
form. For Mode I stress intensity factors at an imbedded 
crack tip it was possible to obtain a two-digit accuracy 
without any difficulty. However, in the calculation of 
the stress intensity factors at the interface crack tips 
and at the intersection of the crack and the interface 
there were convergence difficulties. 
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Figure 5. Stress intensity factor in bonded layers with 
a pressurized edge crack. p-|=U2/3> v-|=V 2 = 0.3. 





t-igure b . The effect of the stiffness ratio y-j/y^ oh the 
stress intensity factor in Bonded layers with 
a pressurized edge crack . v-j * 3 * 

2l2~ b 2 = ^2' 




h , / h 2 

Figure 7. The effect of thickness ratio on the stress 

intensity factor in bonded layers with a pres- 
surized edge crack. y-j = 3u2’ v^=0.3, V2=0.25, 

$2=0. 385339 (for the crack touching the inter- 
face: h«). l~<h 9 : k(b ? ) = lim V2(x 2 -b 2 )a 2 (x 2 ,0) , 

$ x 2^ b 2 
£ 2 = h 2 : k ( b 2 ) = 1 i m /2 x ] 2 a lyy ( x 1 ,0) . 





Figure 8. Same as Figure 7. u-^Ug/S, v-j=0.25, v 2 =0.3 
3 2 =0. 624348 (for ^ 2 =h 2 ). 


k(efe) 


Stress intensity factors for an imbedded crack 
touching the interface. y-| = 3^ » =0 . 25 , ^2=0. 3, 

external load: uniform crack surface pressure. 
k(a 2 ) = l im o ^ y y ( x 2 » ^ » k ( b 2 ) — 1 x -j' ■ o 1 ^ 



Figure 10. 


Stress intensity factor in bonded layers con- 
taining an edge crack and subjected to uniform 
bending away from the crack region . v i ' 3 * 
solid lines: u-j = 3u 2 > dashed lines: 
k(b 2 )-lim /2l x2~t>2 )a 2 yy (x 2 »0) . 




h | / h 2 

Figure'll.. Stress intensity factor in the composite plate 
with a broken layer subjected to uniform bend- 
ing. v 1 =v 2 =0.3. Power of stress singularity 
6 2 =0. 400470 for p-j=3y 2 and B 2 =0. 620492 for 
u-i=Uo/3. k(b 0 ) = lim /2 x-, ° i x/ «, ( x i » ^ ) • 
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Figure 12. Stress intensity factors in bonded layers with 
double edge cracks . u-| = 3 ua * v-j =V 2 = 0 . 3 , exter- 

nal loads: uniform crack surface pressures, 
p -j = 3p 2 > h-|=h 2 = h, 21 ^ 21 ^ 21 . k ( a -j ) . * 
i im >/2 ( a, -Xt )a, .... ( x-, ,0 ) , k (b 0 ) = lim/2(x 0 -b 0 )a, 



Figure 13. 


Stress intensity factors for double edge cracks 
Ui “Vo/ 3 » v, =v 9 =0 .3 , h 1 =h 2 , 2£ 1 = h 1 , p 1 = p 2 /3. 




Figure 14. Effect of Poisson's ratio ori the stress inten- 
sity factors ih mm mm With m ddp track 
touching the interface (Hote: a is also depen- 
dent on v, and is give* in tabid l) * 
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Stress intensity factors for an imbedded crack 
crossing the interface. External loads: uni- 
form crack surface pressures p-, , P 2 , p 1 = 3p 2 , 
v-|=V 2 = 0*3, p^ = 3p 2 , h^=h2> b-|=2f.^=h^, k(b-j)- 
lim /2(x 1 -b 1 T a lyy ( x l >°) .k(a2^ = 1in1 ^ a 2' x 2 ^ a 2yy^ x 2. 



a 

Figure 16. Same as Figure 15. k v ~ 1 i m y cr vv (0,y), k * 

a X y+0 XX xy 

1 im y cr (0,y), a=0 .079898. 

y-o xy 


Figure 17. 


Stress intensity factor for an edge crack 
crossing the interface. k(b-,) = 

Urn /5TvT^«, y y(V 0 >- External l0ad5: 
uniform crack surface pressures, P-, . P 2 > 

(l-v?)p 1 /E,.(l-vl)p 2 /E 2 , f/h 2 -1.2, solid 
curve: P 1 =3p 2 , v,=0.3, Vj-0.25, dashed curve 
u,-u,/3, v,-0.25 , \* 2 = 0 . 3 . 



Figure 18. Stress intensity factors at the interface for 

an edge crack crossing the interface, k x = 

1 im y a c (0 ,y ) , k =1 im y^d^uCO >y )> ^/h^=l . 

y->0 xx xy y^-0 y 

u, = 3ii„ % v,=0;3, v„ = 0.2B, a=0. 0601 77, 





h,/h 2 

Figure 19, Same as Figure 18, y-^y^/S, =0 . 25 , V 2 = 0.3, 

a=0. 060177. 




Figure 20. Crack surface displacement for an edge crack 
crossing the interface. Constant crack sur- 
face pressures p-| , P 2 > ( 1 -v\ )P-j /E] - (1 -V 1 ) P 2 / E 2 ’ 
v-| =0 . 35 , v 2 = 0.3, h 1 = 3h 2 , £ ] =0.3h 2 , t ] = (2 x ] /b ] ) - 1 
1 2 — 1 — (x 2 /h 2 ) • 


Stress intensity factors at the intersection 

btr6S . interface in two bonded 

of the free end and the interface n = 

semi-infinite layers. * °xx< 0 ’ ’ W 

lim y“o xy (0.y). u, = 3u 2 . v^S. 3. 02*0 . 25 , a 

.. v n . u . o ~ u o 


, « /* a i n "7 





— = 2 
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Figure 22 


Relative displacement in y-di recti on at the 
free end of two bonded semi -inf ini te layers. 
y,=3y ? , v,=(K3, v 2 = 0.25, t 1 = ( /h 1 ) - 1 , -t| a . 
i*(xJh 2 )i V'M 0 * 01 


Figure 23. Stress intensity factors and the strain energy 

release rate, W for a uniformly pressurized T- 

shaped crack. P 3 =P 2 > P^ 0 ’ h 1 =h 2 =h, b 3 =£ 3 =h, 

y,=3yo» v 1 =Vp=0.3, k,+ik 2 =lim (b 3 -y ) 03 (b 3 +y) a3 
1 L 1 c y-*b 3 

C ff i x x( 0 ’y) +icr ixv (0 ’ y)] ’ w= (k! + k!)/pUi. 





Figure 24. Stress intensity f acton k ( a 2 ) for an internally 
pressurized T-shaped crack (data same as in 
Figure 23) . 
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Figure 25. Stress intensity factors and the strain energy 
release rate W for a uniformly pressurized 
T-shaped crack - the case of broken layer. 
p 3 =p 2 , P 4 =0, u 1 = 3p 2 , v 1 =v 2 = 0.3, ' 




Figure 27. Stress intensity factors for the interface 

crack in bonded layers containing a T-shaped 

crack. p 3 =p 4 =0, u*| “V*2^ 3 ’ v 1 =v 2 = 0 ' 3, £ 3 =h l =h 
(see insert in Figure 26). 
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Figure 29 


Stress intensity factors at the interface 
crack tip in bonded layers containing a T 
shaped crack - the case of broken layer. 


V. 


■V 


, = 0.3, £,-h 9 . 
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Figure 30. 


I 2 

^3/ h 2 

Stress intensity factors k(a 2 ) and k (b 1 ) for 
a symmetric cross-shaped crack in bonded lay- 
ers 1 oaded in tension away from the crack 

region. u-j^Syg* ^ 1 *^ 2 *'®-’ 3» p-j 5 3P|, P 3 = P 4 = 0 » 

h -j =* h 2 ^ 





Figure 31. Stress intensity factors for the interface 

crack in bonded layers containing a symmetric 
cross-shaped crack (data same as in Figure 30). 



k 2 



Figure 32. Stress intensity factors for the interface 

crack in bonded layers containing an internally 
pressurized cross-shaped crack - the case of 


broken layer, y -j ~ 3 y 2 » v i =v 2 = ^ • 3 , • t-j = hi /2 , 
h i - h 2 •» P ^ ~ P 2 * P i ~ 3 • 8 5 7 p £ > P ~ 3 • 





Figure 34. Stress intensity factor k(b-|) in bonded layers 
containing a cross-shaped crack with a broken 
layer and loaded in tension in y direction 
away from the creek region, h-j = h 2 = 2£/3 , v^- 
v 2 =0.3, p 3 = p 4 =0, solid line; u 1 = 3u 2 , p ^ = 3 u 2 , 
dashed line: y, =y 0 /3 , p, =p 0 /3. 
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Figure 35. Stress intensity factors in a bonded semi- 

infinite plate with an edge crack on the inter 

face, l-i i 3 vi 2 ’ ^ i — ^ 2 .~ ^ ^ ^ 3 — ^ 2 5 ^ 2 ^ 2 — ~ ^ 1 ^ 1 5 

M=2p 1 h 1 (h,+h 0 ) . 









